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Isospin is included in the description of Coulomb excitation of multiple giant isovector dipole 
resonances. In the excitation of even-even nuclei, a relevant portion of the excitation strength 
is shown to be associated with 1 + two-phonon states, which tends to be hindered or completely 
supressed in calculations in which the isospin degree of freedom is not considered. We find that the 
excitation cross sections is strongly dependent on the ground state isospin. 

I. INTRODUCTION 



q \ Coulomb excitation of two-phonon giant resonances in heavy ion collisions at relativistic energies was predicted by 
D 1 Baur and Bertulani in 1986 |l| and generated considerable interest during the last few years [§d4[. The double giant 
£^ \ isovector dipole resonance (DGDR) has now been observed in several nuclei: 136 Xe H, 197 Au |q|, and 208 Pb 
_ -. • The double giant isoscalar quadrupole resonance has been also observed in the 40 Ca ( 4 %3a, 40 Ca pj reaction, at 44 • A 
MeV laboratory energy Q . Data on the DGDR were confronted with results of coupled channels Coulomb excitation 
y—( , theory |l(],[ll|] and the major conclusion reached was that theoretical predictions underestimate the data by a factor 
' 2-3 in the cases of 136 Xe and 197 Au. A similar discrepancy, albeit appreciably smaller, was found in 208 Pb fll2| . 

Several effects, not taken into account in the coupled channel theory, were considered to explain this discrepancy. 
Among these are anharmonicities |]l 3|J 1 4|| and the Brink- Axel mechanism jl5| , to cite a few. In this paper we examine 
the relevance of isospin effects in the excitation process relying on the pionee ring work of Fallieros, in which the 
isospin splitting of the isovector giant dipole resonance was first analyzed fl6Hlq|. We next present an extension 
of this work to the double giant isovector dipole states |l9| , before turning to more technical details of the coupled 
channels calculation. 

""^ ' The relevant dipole excitation operator is the component T3 = of an isovector (T = 1) operator, Dt=i,t 3 =o- Thus, 
when acting on a target nucleus with isospin quantum numbers T, T 3 = T, two GDR states will be generated having 
1 1 ■ isospin quantum numbers Tj = T and Tj = T + 1 (Tj = T — 1 is forbidden due to charge conservation) , where the 
label I has been introduced to designate the (intermediate) one-phonon state. If the dipole operator is applied again 
to these one-phonon states, final states with isospin Tf — T, T + 1 and T + 2 will be generated. In order to take into 
account the bosonic nature of isovector phonons in these final states, one must keep track of another isospin quantum 
number, namely the total isospin 3 of the two dipole phonon operator, which can take the values 3 = 0, 1, 2. These 
values will constrain the total angular momentum of the two coupled phonons through symmetry requirements. For 
a nucleus with J v = + ground state, the DGDR may have J 71 " = + , 1 + and 2 + , and for the state 1 + one must have 
Q = 1. If isospin is not taken into account, 1 + states, reached from the GDR, will be quenched ^0|, since by 
itself it is antisymmetric under exchange of the two phonons. However, 1 + states will in general contribute to the 
excitation cross section, if explicit reference to its Q = 1 nature is made. In this case, the exchange symmetry is odd 
both in the spin and isospin spaces, so that the product has even symmetry, as required. 

The energy splitting of the isodoublet GDR was studied in |ll| and was found to be related to the symmetry energy 
and to the average particle-hole interaction, leading to the estimate 

A« 1 = sgi - Eft = 60 [MeV] . (1.1) 

The energy splitting of the isotriplet DGDR can be estimated in a similar way. Since the Tf = T + 2 state is a double 
isospin analog state it involves twice the displacement energy, and we may write 



= 4+2 ~ 4 2) = 120 ^ [MeV] . (1.2) 



Furthermore, for the Tt = T + 1 state we may write 



A^L = 4+1 - e t' - 60 1 ~r 1 t MeV ] • ( L3 ) 



T+l 



1 



We give in table I the energies of the isospin doublet and triplet states in 208 Pb (T = 22) as obtained from these 
expressions. 



Energy shift (MeV) 


^u» pb ( T = 2 2) 


4 »Ca (T = 4) 


A (l) 

A (2) 

A (2) 
T+l 


6.63 

13.85 

6.63 


6.25 

15 

6.25 



Table I. Isospin splitting (in MeV) of one- and two-phonon states in 208 Pb and 48 Ca. 



Having given the above account on the isospin structure of the one and two-phonon dipole states, we now turn to the 
required modifications of the coupling interaction for the coupled channels calculation. 



II. EXCITATION OF MULTIPHONON STATES 

A. The coupling interaction 

The coupling interaction for the nuclear excitation i — > / in a semi-classical calculation for an electric (tt = E), or 
magnetic (tt = M), multipolarity, is given by (eqs. (6-7) of ref. pi]]) 

Wc = — = E^v(r), (2.1) 



where 



(r) = ( -l) A+1 |£ 1 J j^^Q,x» ({, r) M (ttX, p) . (2.2) 

In this expression b is the impact parameter, r = ^yvt/b is a dimensionless time variable with 7 = (1 — /3 2 ) and 
(3 = v/c being the usual relativistic parameters. The energy scale is set by eo = jhv/b and the quantity Q-nX^ (£) r )> 
with £ defined as the adiabatic parameter £ — £,if = {Ef — -Ei) /eo? depends exclusively on the properties of the 
projectile-target relative motion. The multipole operators acting on the intrinsic degrees of freedom are, as usual, 

M{E\,n) = J d 3 r p(r) r A ^(r) , (2.3) 

and 

M(Ml,|i) = ~ J d 3 r J(r).L (rY^) , (2.4) 

We treat the excitation problem by the method of Alder and Winther |2^| . We solve a time-dependent Schrodinger 
equation for the intrinsic degrees of freedom in which the time dependence arises from the projectile-target motion, 
approximated by the classical trajectory. For relativistic energies, a straight line trajectory is a good approximation. 
The wave function is expanded in the relevant eigenstates of the nuclear Hamiltonian, {| fc); k = 1, N}, N being the 
number of relevant intrinsic states included in the coupled-channels (CC) problem. 

In order to simplify the expression for the coupled equations we define the dimensionless parameter rjy by the 
relation 



The coupled channels equations can then be written in the form p] 



2 



da k (r) 
dr 



N 
r—1 



In what follows we concentrate on the El excitation mode. In this case, we have 

i 2 

. C/ 



^ 3 (*)-ifQ)%(*) 



Qex±i(£,t) =T0 3 (t) 



(2.6) 



(2.7) 



1 /2 

where </>(r) = (l + t 2 ) 

The excitation probability Pj(b) of an intrinsic state | j) in a collision with impact parameter b is obtained from 
the amplitudes Oj(r) at asymptotically large times, in terms of an average over the initial and a sum over the final 
magnetic quantum numbers. The cross section is then obtained from the classical expression 



o j = 2ir J P 3 {b) T{b) bdb . 
where the impact parameter dependent transmission factor T(b) accounts for absorption pj| 



(2.8) 



B. The 1-phonon states 



Consider the excitation of a nucleus with ground state spin zero (any even-even nucleus) and isospin To, with 
its 3-component T3 = Tq. In terms of the relevant quantum numbers, these states are written as: \j >= 
\Ej n ^;JjMj;TjT3j >, where n is the number of phonons, E the energy, Jj and Mj are respectively the spin and 
its z-component quantum numbers, and Tj and T$j are the isospin and its third component. Note that due to charge 
conservation, all states have T 3j = T . The ground state and the 1-phonon states are given in the Table II below. 



n 


E 




T 








0+ 


To 


1 


Egdr 

EgDR + &Ta+l 


1- 
1- 


To 

To + 1 



Table II. Ground-state and one-phonon states with angular momentum and isospin dependence. 



The energy splitting A^P , which is assumed to depend exclusively on isospin, is given by eq. 

In order to calculate the matrix elements Mki (7rA, [i) between initial, i, and final states, k, we use the Wigner- 
Eckart theorem in spin-isospin space and (except for the energy dependence) assume that the reduced matrix elements 
are isospin independent. We get 



1.1 



A4f (^) = (E£;lM k ;T k T M (El,/j) 



eW;Q0;T o T o 



1 



1 m 1 T 



-M k n 



To Ti 



_ ,_^l-Mf+T k -T 

x{l\\M{El)\\Q) . (2.9) 
The value of the reduced matrix element (1 1 1 A^f (El) ||0) can be extracted from the energy-weighted dipole sum-rule 



S = 12 ( E k ] - Ei0) ) (E%;lM k ;T k T M(El,n) E<®- t 00; T T 



M k T k 



3 h 2 NZ 
4n2mN A 



(2.10) 



Inserting the matrix-elements of eq. 2.9 in eq. 2.10, we get 



5 = E(^ 1) - S(0) ) w\ M ( M )H°) 



2 (T k 1 T 
T T 



AI, 



1 

-M k ^ 



(2.11) 



3 



Using the relation 24 1, 



E 

M k 



1 1 o \ 2 (\ 1 o\ 2 1 



ii n o y 3 1 



(2.12) 



the sum rule takes the form 



5 = < I . To o T ) Egdr + I -To T ) r GI3it + A ^o+i 



(2.13) 



Using the explicit forms of the Wigner coefficients 



To 1 T 
-T T 



To + 1 1 T 
-T T 



we obtain the reduced matrix element 



T n 



"I 1/2 



(T + 1) (2T + 1) 
1 



"I 1/2 



(T + 1) (2T + 3) 



(2.14) 



(1\\M (251)||0> 



35 



EgdrF(T ) 



(2.15) 



with 



F(T ) = 



1 



T n 



T + 1 2T + 1 



A (D 
Egdr 



1 



2T + 3 



(2.16) 



C. The two-phonon states 

The two-phonon states must be symmetric with respect to the exchange of the two phonons in spin and isopin. 
This symmetry can be tracked by using the coupling scheme KI1I2) 2To; TfTo). One has to distinguish the two cases 
3 = (isospin even) and 3 = 1 (isospin odd). The states which are spin-isospin symmetric correspond then to the 
combinations 



(2^30 = (0,0); (0,2); (2,0); (2,2); (1,1). 



(2.17) 



(2) 

The two-phonon states are represented as \j >— |25) ; JjMj] (I1I2) STo 2/ To >. We list the main quantum numbers 
of the two-phonon states in table III. 
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To 
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To 
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T - 


- 1 
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Edgdr 




0+ 


T 






Edgdr 


hA (2) 


2+ 


To 






Edgdr - 


0+ 


T - 


hi 




Edgdr - 


- A (2) 


2+ 


T - 


hi 




Edgdr - 


- A (2) 


0+ 


T - 


h2 




Edgdr - 


- A (2) 


2+ 


T - 


-2 



Table III. Two-phonon states with angular momentum and isospin dependence. The isospin dependence arises from 

the coupling of the phonon isospins 3 = li ® 1%. 
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The isospin shifts A' 2 ' in the above table are given by eqs. 1.2 and [i~3| . 

To obtain the excitation amplitudes for the one phonon to two phonon transitions we need to calculate the following 
matrix elements: 



M% ] (El,n) =^W;J / M / ;(li,l 2 )93b;T / r M(El, l i) 



A 2 ). 



E^;J k M k ;T k T 



(2.18) 



This can be done by changing the final state coupling scheme. We use the relation (ref. (g4j],eq. 6.1.3, p. 91) 
\JM; (lil 2 ) 3T ; T f T ) = ]T (-lf l+Tf ^/(2Z + 1) (27] + 1) 

Ti 

}, y | j|JAf;(l 2 (l 1 To)T i );r / r ) . 

J ' J 



(2.19) 



and obtain 



(£1, ») = J2 (-l) Ti+T/ V(23+1)(2T, + 1) 



1 1 3? 



E^-JjMf, (l a (l^o) T,) ; T/T (£1, /x) 



E^;J k M k ;T k T 



(2.20) 



We next use the Wigner-Eckart theorem in spin-isospin space, assuming that the reduced matrix element are spin 
and isospin independent and vanish unless T k = Tj. We get 

(E^JjMfi (1 2 (IiToJTOjT/ToI X (M,^) J fc M fc ; (liT )T fc T ) = 
W • (-1)^-^— (_ M/ J 2) ( T/ To S) (2 ' |A1 ' 



(2.21) 



Using {2\\M (El) ||1) = V2(1\\M (E1)\\0) in cq. [E2l], eq. ^ becomes 

TW^ 15 (.El, /x) = 72(-l) J ^ M ^ +2T ^ To+Tfc (1||7W (Sl)||0) A /(29 + l)(2T fe + l) 
'1 1 0\ /T/ 1 T k \ f 1 1 3 



M f nOl V -T T M T 2) T fc 



(2.22) 



III. APPLICATIONS 



In this section, we apply our results to the excitation of one- and two-phonons states in 208 Pb and Ca target nuclei, 
in collisions with relativistic 208 Pb projectiles. Before we present the numerical results, we rewrite the matrix elements 
so that they can be used as input to the coupled-channel code RELEX for Coulomb excitation plfl . They become 
"effective reduced matrix elements" that incorporate the isospin dependence of -M^P (ttX, fx) (eq. 2.9). Namely, 



(l||A4(i51)||0>&"> 



(-1) 



T fc -To 



while for equation 2.22 it is more convenient to define 

(2\\M (El)\\l)% ff) = V2(-l) 2Ti - To+Tk (l\\M (£1)||0)V(23 + 1) (2T k + 1) ^ 



1 T fe 
-To T 



1 1 Q? 
T T f T k 



(3-1) 



(3.2) 



The modified reduced matrix-elements for one- and two-phonon excitations are presented in tables IV and V, in the 
cases of 208 Pb and 48 Ca, respectively. For comparison the original reduced matrix-elements are also given in the last 
column. 
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27 
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22 
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824 
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33.6 







23 


23 
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-6.56 


20.1 


33.6 




2 


23 


23 


2 


-6.56 


20.1 


40.8 







23 


24 


2 


-2.83 


20.1 


40.8 




2 


23 


24 


2 


-2.83 



Table IV. Reduced matrix elements for transitions from the GS to one-phonon states and also from one to 
two-phonon states, for 208 Pb in units of fm.e. The numbers within parentheses in the last column are the 
corresponding reduced matrix elements without consideration of isospin. To avoid ambiguities in the two-phonon 
final states, we also indicate the isospin of the coupled-phonon pair. 



Ek 


Ef 


Jk 


Jf 


T k 


Tf 


3? 


(n+l||.M (.El)||n) 





19.2 





1 


4 


4 




2.66 (3.00) 





25.4 





1 


4 


5 




1.2 


19.2 


38.4 


1 





4 


4 





-2.17 (4.24) 


19.2 


38.4 


1 


2 


4 


4 





-2.17 (4.24) 


19.2 


38.4 


1 


1 


4 


4 


1 


0.594 


19.2 


44.6 


1 


1 


4 


5 


1 


-1.08 


19.2 


38.4 


1 





4 


4 


2 


3.01 


19.2 


38.4 


1 


2 


4 


4 


2 


3.01 


19.2 


44.6 


1 





4 


5 


2 


1.31 


19.2 


44.6 


1 


2 


4 


5 


2 


1.31 


19.2 


53.4 


1 





4 


6 


2 





19.2 


53.4 


1 


2 


4 


6 


2 





25.4 


38.4 


1 





5 


4 





1.09 


25.4 


38.4 


1 


2 


5 


4 





1.09 


25.4 


38.4 


1 


1 


5 


4 


1 


1.19 


25.4 


44.6 


1 


1 


5 


5 


1 


-2.15 


25.4 


38.4 


1 





5 


4 


2 


0.548 


25.4 


38.4 


1 


2 


5 


4 


2 


0.548 


25.4 


44.6 


1 





5 


5 


2 


-1.76 


25.4 


44.6 


1 


2 


5 


5 


2 


-1.76 


25.4 


53.4 


1 





5 


6 


2 


-1.92 


25.4 


53.4 


1 


2 


5 


6 


2 


-1.92 



W77 



Table V. Same as Table V, for 48 Ca 
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The calculated excitation cross sections for one and two phonon states in the 208 Pb (650 MeV • A) + 208 Pb collision 
are given in table VI. They are also plotted in hgure la, as a function of the excitation energy. In this hgure, the 
dominant spin and parity is indicated in each case. For a comparison, corresponding results neglecting isospin are 
given within parentheses in table VI and are shown in figure 16. If isospin is neglected, only T — 22 (the ground 
state isospin of 208 Pb) states are populated. Namely, the GDR state at 13.5 MeV and + and 2 + DGDR states at 
27.0 MeV. With the inclusion of isospin, about 3 % of the GDR cross section is associated with the population of the 
T = 23 state at 20.1 MeV, as can be seen in table VI and in figure la. This corresponds to the exhaustion of about 
6.3 % of the GDR sum rule. The consequences of the isospin degree of freedom on the DGDR population are more 
important. Although most of the cross section remains associated with the T = 22, J 71 = + and 2 + states at 27.0 
MeV, 7 % of the total DGDR cross section then arises from the population of the T = 23 states at 33.6 MeV, of which 
over 90 % corresponds to the non-natural parity J 71 = 1 + state. Thus, 6 % of the DGDR cross section goes to the 
excitation of the J 71 = 1 + state, which would be forbidden in the usual harmonic oscillator picture (without isospin). 
On the other hand, some calculations using RPA descriptions of the giant resonances find non-vanishing population 
of such states. Lanza et al. |D| find negligible populations while Bertulani et al. get about half of that of the 
present work. This result could be traced back to the fact that their second order transitions cancel exactly, so that 
1 + states can only be reached through higher order coupled channel processes. 

Table VII and figure 2a give similar results for the excitation of GDR and DGDR states in 48 Ca, in the collision 
208 Pb (650 MeV ■ A) + 48 Ca. Corresponding results neglecting isospin are given in the same way as above. For this 
system, isospin plays a more important role due to the lower isospin quantum number (T = 4) of the 48 Ca ground 
state. In this case, the dominant T = 4 GDR state at 19.2 MeV loses more than 10 % of its cross section to the T = 5 
l - state at 25.4 MeV, which exhausts about 21 % of the GDR sum rule. Moreover, the DGDR cross section is very 
much affected by isospin. Figure 2a indicates that the cross section for T = 5 DGDR states at 44.6 MeV reaches 32 % 
of that for the dominant T = 4 DGDR states at 33.6 MeV. It is also important to discuss the J 71 distribution of the 
DGDR cross section. Since over 95 % of the T = 5 DGDR cross section corresponds to J n = 1 + , the population of 
states with this spin and parity is rather large. 



E (MeV) 


J* 


T 


3 


a (mb) 




GDR 13.5 


1- 


22 




2240 


(2597) 


20.1 


1- 


23 




59.86 




DGDR 27.0 


0+ 


22 





4.489 


(15.80) 


27.0 


2^ 


22 





31.16 


(117.5) 


27.0 


1+ 


22 


1 


0.8608 




33.6 


1+ 


23 


1 


8.373 




27.0 


0+ 


22 


2 


10.85 




27.0 


2 + 


22 


2 


69.89 




33.6 


0+ 


23 


2 


0.1257 




33.6 


2 + 


23 


2 


0.3732 




40.8 


+ 


24 


2 


0.0237 




40.8 


2+ 


24 


2 


0.0874 





Table VI. Excitation cross sections of one and two phonon states in the collision 208 Pb (640 ■ A MeV) + 208 Pb. 



E (MeV) 


J w 


T 


3 


a (mb) 




GDR 19.2 


1- 


4 




318.5 


(405.5) 


25.4 


1- 


5 




33.11 




DGDR 38.4 


0+ 


4 





0.2842 


(2.137) 


38.4 


2+ 


4 





0.4453 


(3.613) 


38.4 


1+ 


4 
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0.6695 




44.6 


1+ 


5 


1 


1.275 




38.4 


0+ 


4 


2 


1.117 




38.4 


2+ 


4 


2 


1.640 




44.6 


0+ 


5 


2 


0.0155 




44.6 


2+ 


5 


2 


0.0366 




53.4 


0+ 


6 


2 


0.0268 




53.4 


2+ 


6 


2 


0.0375 





Table VII. Same as Table VI, for the excitation of 48 Ca in the collision 208 Pb (640 ■ A MeV) + 4s Ca. 
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IV. CONCLUSIONS 



In this paper isospin is taken into account in the excitation of the double giant dipole resonance. We have used 
a semiclassical coupled-channels formalism to calculate excitation probabilities and cross sections for the collisions 
208 Pb (640 • A MeV) + 208 Pb and 208 Pb (640 • A MeV) + 48 Ca. We have assumed that the eletromagnetic matrix 
elements are isospin independent and adopted the isospin splitting of energy levels as given by Akyiiz and Fallieros 
||l8f . It has been shown that isospin leads to a redistribution of the strengths of the electromagnetic matrix elements 
such that the probability for Coulomb excitation of J v = 1 + DGDR states is enhanced. This enhancement depends 
on 3J and 6J coefficients in isospin space so that it becomes more relevant for nuclei with low ground state isospin. 
Consequently, this effect turned out to be much stronger in the excitation of 48 Ca than in the excitation of 208 Pb. 
One should then expect that isospin splitting should contribute to make de DGDR broader, particularly for nuclei 
with low ground state isospin. This result suggests that our formalism should be especially relevant to study the 
excitation cross sections of a family of isotopes with charge numbers varying from the proton to the neutron dripline. 
Study along these directions is underway. 
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